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1 Introduction 

In many extensions of the standard model (SM) of particle physics [1-3] the spontaneous break- 
ing of the electroweak (EW) symmetry |4]-[9l is associated with new strong dynamics appear- 
ing at the TeV scale. For instance, the origin of the new dynamics may be due to new inter- 
actions flT0"llT2"J , compact extra dimensions f|T3l fl4 | . or a composite Higgs boson [TBI [16|. In 
such scenarios the SM is an effective low-energy theory valid for energies smaller than a new- 
physics scale A. In these theories, one expects the existence of new resonances coupling to pairs 
of vector bosons (ZZ, WZ, and WW). A minimal ultraviolet completion of this effective theory 
for composite models is described in Ref. [17]. Other examples include Randall-Sundrum (RS) 
gravitons Gkk H3lH4l coupled to ZZ and WW, or technimesons IIT81IT9I coupled to WZ. Limits 
from previous searches and from indirect bounds (e.g. in the EW sector and flavor physics) 
place the masses of these proposed RS resonances at or above the TeV scale fl20Tl26l . These 
scenarios could be tested at the Large Hadron Collider (LHC), as long as A ~ O(TeV), as sug- 
gested by the EW symmetry breaking scale. This analysis is sensitive to searches for resonances 
starting at 700 GeV and above. However, there are other theories that predict light resonances 
(e.g. low-scale technicolor) |fl8l[l9| . 

In this Letter we present a search for heavy resonances decaying to WZ and ZZ final states, with 
one boson being a Z decaying to leptons, namely Z — > £ + £~ {I = }i, e) or vv, and the second bo- 
son decays to hadrons, V (V=W, Z) — > qq. For heavy resonances the decay of each V produces 
a highly boosted system in which the two fermions are emitted within a small opening angle 
in the laboratory frame. The hadronization of the V— >■ qq quarks would then produce two par- 
tially overlapping jets reconstructed as a single jet with mass close to the V mass, a topology 
very different from that of a typical quark or gluon jet. Monte Carlo (MC) simulations suggest 
that more than ~70% of the decays would produce a merged-jet topology for resonances heav- 
ier than ~800 GeV. This feature is exploited in a VZ final state, to discriminate a possible signal 
from the SM background (mainly coming from Z+jets events). 

Thus, in this study we consider three final states: one heavy jet and either Z — > e + e~, Z — > 
, or Z — > E™ ss , where E™ 1SS is the characteristic signature of neutrino production. We 
characterize the signal as a peak in the invariant mass of the VZ system (transverse mass in 
the case of Z — > vv decays). The search is performed with a data sample of pp collisions 
corresponding to an integrated luminosity of 5.0 fb _1 collected by the Compact Muon Solenoid 
(CMS) detector at the LHC at y/s = 7 TeV in 201 1 . 

Results are presented in terms of two benchmark scenarios: i) the Sequential Standard Model 
(SSM) in which a new gauge boson W with the same couplings as the SM W boson decays to 
a WZ pair; ii) a RS graviton, Gkk/ decaying to ZZ. In both scenarios we search for resonances 
heavier than 700 GeV, where the considered boosted topology becomes relevant. The ratio of 
the 5-dimensional curvature to the reduced Planck mass (k/ Mpi), which acts as the coupling 
constant in the RS model, is typically used as the phenomenological parameter in RS graviton 
searches. For the RS graviton scenario we consider values of the coupling parameter k/ Mpi up 
to 0.3. 

Previous searches have been carried out in the context of both the SSM W' and RS graviton 
theoretical models. The most stringent limits have been produced at the LHC by the ATLAS 
and CMS collaborations in a large number of final states: W' tv EZHHI, W' -> tb E9ll30l , 

W' ->• WZ -> 3£v EIH32I, and Gkk -> U Q21I32I, Gkk -> 77 HUM and Gkk -> ZZ ->■ Ujj 

Eg [an. 
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3 Collision and Monte Carlo Datasets 



2 The CMS Detector 

Here a brief description of the CMS detector is given with an emphasis on the elements most 
relevant for this analysis. A more detailed description can be found elsewhere [37|. A cylin- 
drical coordinate system about the beam axis is used, in which the polar angle 9 is measured 
with respect to the counterclockwise beam direction and the azimuthal angle, (p, is measured 
in the x-y plane in radians, where the x axis points towards the center of the LHC ring. The 
quantity r\ is the pseudorapidity defined as t] = — ln[tan(0/2)]. The layout comprises a super- 
conducting solenoid providing a uniform magnetic field of 3.8 T. The bore of the solenoid is 
instrumented with various particle detection systems. The inner tracking system is composed 
of a pixel detector with three barrel layers at radii between 4.4 and 10.2 cm and a silicon strip 
tracker with 10 barrel detection layers extending outwards to a radius of 1.1 m. Each system is 
completed by two end-caps, extending the acceptance up to \rj\ < 2.5. A lead-tungstate crystal 
electromagnetic calorimeter with fine transverse (A//, Acp) granularity and a brass /scintillator 
hadronic calorimeter surround the tracking volume and cover the region \ tj\ < 3. CMS also has 
extensive forward calorimetry. The steel return yoke outside the solenoid is instrumented with 
gas-ionization detectors which are used to identify muons in the range \r}\ < 2.4. The barrel 
region is covered by drift tube chambers and the end cap region by cathode strip chambers, 
each complemented by resistive plate chambers. 

3 Collision and Monte Carlo Datasets 

The preselection of the datasets for the analysis is different for the "dilepton" (VZ— > qq £ + £~, £ = 
e, }i) and the "E mi8S " (VZ— > qq vv) channels. For the dilepton channels, we consider events that 
were recorded with double-electron or single-muon triggers. The trigger thresholds changed 
with time, as a consequence of the increasing peak luminosity and the changes in running 
conditions. The tightest thresholds used in the trigger (i.e. 40GeV for the single-muon trig- 
ger and 17GeV for the dielectron trigger) are looser than the corresponding offline analysis 
requirements. Typical trigger efficiencies exceed 83% (95%) for the electron (muon) triggers. 
For the E™ 188 channel, we use triggers requiring at least one calorimetric jet and missing trans- 
verse energy. These triggers have efficiencies of more than 99% for events with a leading jet of 
transverse momentum pj > 160 GeV and E™ ss > 300 GeV after offline reconstruction and cor- 
rections, which allows resonances heavier than 1000 GeV to be probed with an efficiency above 
20%. We use MC samples to study the signal and background. We consider the SM back- 
ground processes that could contribute with two leptons and a (massive) jet in the final state. 
The summary of the signal samples is given in Table [lj and the background samples, in Table [2] 
The pythia 6.424 [38] leading-order (LO) generator with tune Z2 fl39| is used to generate the 
signal events and simulate the parton showering, with a full simulation of the detector based 
on Geant4 9.4 package. Mass-dependent K factors are applied. For the Gkk analysis, next- 
to-leading order (NLO) corrections are calculated using the "two cutoff phase space slicing" 
method ll4"0"ll4"T| in the diphoton final state. For the W analysis, the next- to-next-to-leading or- 
der (NNLO) corrections are calculated with fewz [42j in the leptonic final state. These K factors 
are used for lack of better (N)NLO calculations for the final states considered. The background 
samples are generated using the MadGraph 5.1.1.0 matrix-element generator Il4"3"ll44"| , while 
PYTHIA is used for the parton showering and hadronization, with the same version and tuning 
as for signal samples. The parton distribution function (PDF) used is CTEQ6L1 [45J. Jets are 
matched to partons using the MLM scheme [46 J. 
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Table 1: Signal Monte Carlo samples. The listed cross sections are pythia LO, per channel 
(e + e~ or or vv). The notation vv includes all three neutrino flavors. The K factors com- 

prise NLO (NNLO) corrections for the Gkk (W ) samples. The Gkk samples are generated with 
k/M n = 0.05. 



Mass (GeV) 


Cross Section (pb) 


K factor 






-> qqfi^i (e^e or ji 


V- ) 


Gkk -> qqvv 




750 




O.OD X 1U 




D.VO X 1U 


1.75 


1000 




l.DZ X 1U 




q no v- 1 n — 3 


1.78 


1250 




3.47 xl0~ 4 




2.16 xlO" 3 


1.79 


1500 




8.83 xl0~ 5 




5.24xl0- 4 


1.78 


1750 




3.43 xl0~ 5 




2.04 xlO" 4 


1.76 


2000 




7.05 xl0~ 6 




4.18 xl0~ 5 


1.76 




w'- 


> qql + £~ (e + e~ or ^ 


r) 


W qqvv 




700 




1.267xHT 2 




7.45 xl0~ 2 


1.35 


800 




6.815xKT 3 




4.06xl0- 2 


1.35 


900 




3.842 xKT 3 




2.31 xlO" 2 


1.34 


1000 




2.286 xKT 3 




1.39 xl0~ 2 


1.33 


1100 




1.413xKT 3 




8.45 xlO" 3 


1.32 


1200 




8.955 xKT 4 




5.34 xlO" 3 


1.31 


1300 




5.750X10- 4 




3.42 xlO" 3 


1.30 


1400 




3.784 xKT 4 




2.27xl0" 3 


1.28 


1500 




2.554 xKT 4 




1.50 xl0~ 3 


1.26 



4 Reconstruction and Event Selection 

Events are required to have at least one primary vertex of good quality, where the vertex is re- 
constructed within ±24 cm of the nominal interaction point along the beam axis, with a trans- 
verse distance from the beam spot of less than 2 cm fl4"7|. The events are reconstructed with the 
particle-flow (PF) technique [48J . The PF algorithm reconstructs a complete list of particle can- 
didates in each event from the measurements in all the components of the CMS detector in an 
integrated fashion. The algorithm separately identifies muons, electrons, photons, charged and 
neutral hadrons. Charged hadrons that are consistent with primary vertices other than the lead- 
ing one (defined as the vertex with the largest sum of track pj) are removed from the collection 
of particle candidates used to reconstruct the jets, to mitigate the effects of multiple proton- 
proton interactions within the same bunch crossing (pileup). Electrons are reconstructed as 
isolated objects in the calorimeters which satisfy requirements on the shower shape and the ra- 
tio of the hadronic to the electromagnetic energy deposits. Due to the boosted topology of this 
analysis, some care is needed when reconstructing the Z — > e + e decay: each reconstructed 
electron interferes with the isolation definition of the other electron and has to be excluded from 
the isolation calculation in order to avoid introducing inefficiencies. The isolation criterion for 
electrons is the combined relative isolation Rj so built upon information from the tracker, ECAL 
and HCAL. In calculating Ri S0 , the track momenta and energy deposits, excluding those asso- 
ciated with the electron itself, are summed in a cone of radius AK < 0.3 around the electron 
direction, where AR = (A(p) 2 + (At]) 2 , and divided by the electron transverse momentum. 
Muon tracks are built by combining a track from the inner tracker and a track from the outer 
muon system. No explicit isolation requirement is imposed on the muon candidates. Lepton 
(electron and muon) candidates are required to have a transverse (longitudinal) distance to the 
leading vertex smaller than 2 (5) mm. Jets are clustered from the reconstructed PF particles us- 
ing the infrared-safe anti-fcx [49] algorithm with a distance parameter of 0.7, as implemented in 
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4 Reconstruction and Event Selection 



Table 2: Background Monte Carlo samples. The notation I stands for electrons, muons, or taus. 
The notation vv includes all three neutrino flavors. 



Channel 


Cross Section 
(pb) 


Simulation Details 


Dilepton Channels 


W+jets 

it 

7V+jets 

Z/7*(£+£-)+jets 
WO)W(»+jets 
W(qq)Z(£+£-)+jets 
Z(qq)Z(^+£-)+jets 


212.5 
157.5 
56.5 
25.1 
3.8 
1.14 
0.57 


LO(p™ > 100 GeV) 
NLO 
LO 

LO(pf > 100 GeV) 
LO 
LO 
LO 


E™ 188 Channel 


QCD multijets 
QCD multijets 
Z/y(vv) + jets 


5856.0 
122.6 
32.92 


LO (500 <H T < 1000 GeV) 
LO (H T > 1000 GeV) 
LO (Hj > 200 GeV) 



FAST JET 11501 |51~H . The jet momentum is determined as the vector sum of all particle momenta 
in the jet, and is found in the simulated data to be within 5% to 10% of the true momentum of 
the generator-level jet over the whole pj spectrum and detector acceptance Il52l . An area-based 
correction is applied to take into account the extra energy clustered in jets due to additional 
proton-proton interactions within the same bunch crossing, and for the average effect of out- 
of-time pileup interactions |53ll54| . Jet energy corrections are derived from the energy balance 
of dijet and photon+jet events. Additional selection criteria are applied to each event to re- 
move spurious jet-like features originating from isolated noise patterns from the hadron or the 
electromagnetic calorimeters. The offline missing-transverse-momentum vector (p~r miss ) is cal- 
culated as the negative vector sum of the transverse momenta of all PF particles reconstructed 
in the event, and its magnitude is denoted by E™ 188 . 

4.1 Dilepton channels 

Candidate events are required to have at least two good quality reconstructed leptons within 
the detector acceptance (\rj\ < 2.5 for electrons and \t]\ < 2.4 for muons, with at least one muon 
within \tj | < 2.1 at the trigger level) with pj > 45 GeV. We also require at least one jet in the 
event reconstructed with pj > 30 GeV within \tj\ < 2.4. 

Whenever two same-flavor, opposite-sign leptons are found in the event, a Z candidate is 
formed summing the four-momenta of the leptons. We select the Z candidates by requiring 
their invariant mass to be in the range 70 < Mz < 110 GeV and with a transverse momentum 
Pj > 150 GeV. If there are multiple Z candidates, the one with mass closest to the nominal Z 
mass is selected. The requirement that the dimuon mass is consistent with a Z — > }i + p.~ decay 
strongly suppresses non-prompt muons from jets. 

The V candidate is selected by requiring a reconstructed jet with pj > 250 GeV and \tj\ < 2.4, 
having an invariant mass Mj (computed from the jet energy and momentum calculated as the 
vector sum of the four-momenta of the constituent PF particles) such that 65 < Mj < 120 GeV. 
We require the jet to be well separated from the two leptons forming the Z candidate: AR = 
a/ (A^) 2 + (A(p) 2 > 1.0 for each lepton, where At] (A(p) is the pseudorapidity (azimuthal) dis- 
tance between the jet and the lepton directions. The selection has been optimized by maxi- 
mizing the quantity Ns/^Ns + N-^ (where N$ and Nb are the number of expected signal and 
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E™ ,ss channel 



background events) for the lowest W mass point (700 GeV) considered in this search. 

Once the Z — > ££ and (mono-jet) V— > qq candidates have been reconstructed, we combine 
their four momenta to compute the mass of the parent particle, Mvz- This variable is used to 
evaluate the hypothesis of the signal presence in the data sets analyzed. 

4.2 E™ iss channel 

For the E™ ss channel, background from W-boson decays is reduced through rejection of events 
with isolated electrons or muons with pj > 20 GeV. In order to further reduce leptonic back- 
grounds, we veto on the presence of isolated tracks. For all tracks with pj > 10 GeV and 
\t] | < 2.4, a hollow cone of 0.02 < AR < 0.30 is constructed. The isolation parameter of each 
track is defined as the scalar sum of all tracks with pj > 1 GeV inside the cone, divided by 
the pj of the track. Events containing a track with its isolation parameter smaller than 0.1 
are discarded. Events are then selected if the jet with the highest transverse momentum has 
pj > 300 GeV and \rj\ < 2.4, and E™ 1SS is larger than 300 GeV. In order to reduce the num- 
ber of QCD multijet background events in the signal region, events with more than two jets 
with pj > 30 GeV and \rj\ < 2.4 are discarded. Events with exactly two jets above 30 GeV 
are retained, if the azimuthal angle A(p between the two jets is smaller than 2.8 radians. This 
condition improves the signal over background ratio by reducing the number of QCD dijet 
background events. 

The signal sample is defined as the set of events that meet two extra requirements: the invariant 
mass of the leading jet, Mj, is larger than 70 GeV, and the jet-E™ lss transverse mass, defined as 



M T = y 2 pf Ef ss |l - cos A<£(jet, p™* 

is larger than 900 GeV. Figure TJ on the left (right) shows the two-dimensional Mj vs Mj 
distribution for the simulated SM backgrounds (for a simulated signal sample with Mq kk = 
1250 GeV). 

In contrast to the approach used for the dilepton channels, here we perform a single "event 
counting" experiment by comparing the number of expected background and observed events 
integrated over the region Mj > 70 GeV and Mj > 900 GeV. 

5 Background Estimation 

We are discussing the background estimation separately for the dilepton and E™ ss channels. 
5.1 Dilepton channels 

The analysis of the simulated data shows that the dominant (~90%) background after all selec- 
tion requirements is the inclusive Z production ("Z+jets"), with additional contributions from 
tt+jets and the continuum SM diboson production (WZ and ZZ). The shape and the over- 
all normalization of the expected background Myz distributions are derived from data, with 
additional cross-checks carried out with the inclusive simulated background samples. Effects 
caused by pileup are modeled by adding to the generated events multiple proton-proton inter- 
actions with a multiplicity distribution matched to the luminosity profile of the collision data. 

The background is modeled using a control region consisting of a sideband in Mj (30 < My < 
65 GeV). The remaining selections are applied unmodified to these events, providing a sample 
that is kinematically equivalent to the nominal selection. The robustness of this method against 
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5 Background Estimation 




Leading Jet Mass (GeV) Leading Jet Mass (GeV) 



Figure 1: Distributions of (left) leading jet mass vs. leading jet plus E™ ss transverse mass for 
simulated standard model background sample and (right) RS graviton signal with Mq kk = 
1250 GeV and k/M n = 0.05 sample. 

pileup effects, jet energy scale uncertainties, and variations in the sideband range has been 
confirmed with dedicated studies (Section [6}. 

The procedure is as follows: we first produce the Myz distribution for the sideband selection. 
We define the ratio a (Myz) as the total number of Monte Carlo background entries in the Myz 
spectrum with the nominal (65 < Mj < 120 GeV) and sideband (30 < Mj < 65 GeV) selections: 



N NS (M VZ ) 

a(Mvz) = w 

where Nns(Mvz) is the number of events in the signal region and Nsb(Mvz) is the number of 
events in the sideband region, contained in a bin of the VZ mass distribution centered at a given 
value Myz- We then use the product of the Myz distribution made with the sideband selection 
in the data and the ratio a (Myz) to derive an estimate of the background Myz distribution 
with the nominal selection. Following the example of other resonance searches |55| , we fit the 
a-corrected sideband data Myz distribution to the following analytic function /a (Myz), and 
the fit result is used to parametrize the expected SM background distribution: 



[1- (Ms)]Pi 

/a(Mvz) = P0 r / S , . Mvz . 

^M^+P3log(^P) 

where y/s is the collision energy, p„ i = 0, ...,3 are free parameters of the fit, and Myz is ex- 
pressed in GeV. The fit determines both the shape and the overall normalization of the ex- 
pected background as a function of Myz- The fitting functions are then used to describe the 
expected background in any subregion of the Myz spectrum in the electron and muon chan- 
nels. There are several advantages in using the ratio a (Myz) for the background modeling of 
the Myz distributions: the background estimation becomes insensitive to effects such pile-up 
corrections and integrated luminosity uncertainty which cancel out in the ratio; a (Myz) is less 
sensitive to improper modeling of the matrix element calculation for the background and to 



5.2 E™ ISS channel 
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theory systematics (e.g. normalization and factorization scale, PDFs, etc.) since the background 
composition is similar in the two regions. 



CMS ^s = 7TeV CMS /s = 7TeV 




Figure 2: The comparison of the estimated background (black curve) with the total MC back- 
ground (blue histogram) and the data (black points) for Myz distributions for the electron (left) 
and the muon (right) channels. 

The comparison of the estimated background with the prediction from the simulation and the 
data Myz distributions is shown in Fig. [2j No significant excess of events is observed, with the 
largest deviation appearing in the ~900 GeV region in the muon channel. The tail of the Myz 
distribution, which is the region of interest for the new resonance search, is well described by 
the fit. A discrepancy is observed at low Myz values. Any modeling imperfections, quantified 
as the difference between the best-fit function and the MC simulation estimation, are taken into 
account in the limit calculation by assigning a systematic uncertainty. 

5.2 £!™ ss channel 

By analyzing simulated data we determine that the dominant backgrounds in this channel 
after all selections are inclusive Z — > vv (~70%) and W — > iv (~30%) production, with the 
charged lepton remaining undetected in the latter. To estimate the SM background, we use a 
sideband-based technique similar to that described above, which utilizes events that meet all 
other requirements but the Mj and Mj thresholds. In particular, the events which meet all the 
selection requirements are classified into four regions according to two thresholds in jet mass 
and two thresholds in Mj: 

• Signal region A: Mj > 70 GeV, M T > 900 GeV; 

• Sideband region B: 20 < Mj < 70 GeV, M T > 900 GeV; 

• Sideband region C: 20 < My < 70 GeV, 700 < M T < 900 GeV; 

• Sideband region D: Mj > 70 GeV, 700 < M T < 900 GeV; 

The numbers of events observed in the above regions are denoted as Na, No- 
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5 Background Estimation 



The estimated total background B est in Region A is given by the expression 

N B 1 



Best = Nd 



N c p 



(1) 



where p is a correction factor to account for the correlation between the jet mass and the jet-E™ ss 
transverse mass. The p parameter is estimated from the simulated SM samples by rearranging 
Eq. (TTb in the following way: 

„ ■ Nb 

and setting the values of Na, • • • , Nd to the ones from the SM prediction. Using the values re- 
ported in Table[3]gives p = 0.42 ± 0.02. The value of p thus derived in then reinserted in Eq. (JTJ. 
Setting Nb, No Nd to the yields observed in the data, we obtain an estimate of the remaining 
background B est = 153 ± 20 events. Figures [3] (left) and|3](right) show the comparison between 
the simulated SM background in Region A (scaled to the estimated value B est — a scale factor 
of 11%) and data, together with an example signal for the My and Mj distributions. There is 
agreement between the expected background and data distributions. 

Table 3: Event yields for simulated SM samples, data, and the data/ simulation ratio in the four 
regions described in the text. The quoted uncertainties include those due to the finite statistics 
of the simulated samples. 



Region 


Yield 


Data 


SM Simulation 


Data /Sim Ratio 


A: signal 


N A 


138 


131 ±3 


1.05 ± 0.02 


B: sideband 


N b 


125 


125 ±3 


1.00 ± 0.03 


C: sideband 


N c 


542 


579 ±7 


0.94 ± 0.01 


D: sideband 


N D 


283 


259 ±5 


1.09 ± 0.02 



CMS ^s = 7TeV CMS /s = 7TeV 




Leading Jet Mass (GeV) Jet-E™ ss Transverse Mass (GeV) 



Figure 3: Comparison between p-corrected simulated backgrounds and data in Region A for 
the leading jet mass (left) and jet-E™ lss transverse mass (right) distributions. 
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6 Systematic Uncertainties 

The systematic uncertainties that are considered in this analysis can be divided into two main 
categories: the uncertainty in the determination of the SM background and the uncertainty in 
the expected yields of signal events. All the systematic uncertainties are summarized in Tables|4] 
and for the dilepton and E™ 1SS channels, respectively. The total systematic uncertainty is the 
combination of the signal and background systematic effects, assuming they are completely 
uncorrelated. 

6.1 Background systematic effects 

As we employ a method based on control samples in data for the background determination, 
several systematic effects are eliminated. In the following, we consider the remaining relevant 
uncertainties in detail for the dilepton and E™ ss channels. 

6.1.1 Dilepton channels 

The expected number of background events in each mass window is determined by the integral 
of the function /a(Mvz) i n the corresponding region. The statistical uncertainty is calculated 
by employing the full covariance error matrix of the fit parameter uncertainties in the integral 
of the fitting function in the mass window. The pileup and jet energy scale (JES) systematics 
can potentially affect the background determination through the a(Mvz) ratio and are consid- 
ered separately. The former is found to have a negligible effect. For the latter, the uncertainty 
is evaluated by varying the jet pj to pj ± OjEsiPTrV)/ where <7j ES (p T , t]) is the total jet uncer- 
tainty, and applying the full fitting procedure. The yield differences, in each mass window, 
between the expected background with the positive (N + Bgd) and negative (N_B g d) jet energy 
scale variation with respect to the nominal selection and fit are taken as the ±1 a estimates 
for the JES systematic uncertainty. We also consider several variations in the fitting procedure 
(fitting range, functional form, and sideband definition). These variations are compared to the 
difference in the number of expected background events in the mass window as estimated from 
data and with MC simulation. The largest of the two is used as the systematic uncertainty in 
the background determination. 

6.1.2 E™ iss channel 

To evaluate the robustness of the evaluation of the expected background two tests are con- 
ducted. 

The first test studies the dependence of the correction factor p on the definition of the sideband 
regions. We vary the definition of the sideband regions by changing the thresholds in the My 
and Mj variables in the intervals 20-70 GeV and 650-750 GeV, respectively. We find that the 
resultant variation in the mean estimated background is typically 5% or less, confirming the 
robustness of the sideband method. 

A second test is used to check the propagation of all the uncertainties involved in the B est cal- 
culations. We generate a series of pseudo-experiments with the number of events constrained 
to be equal to that of the actual experiment. We obtain a value of p and calculate the mean 
estimated background in each case. The distribution of the values of B est thus obtained has a 
variance of 20 events. This result is in agreement with the estimated uncertainty on B est ob- 
tained in Section 5.2, using the yields of Nb, Nq and Nd observed in the data. 

The mean expected SM background in Region A, within the uncertainties calculated above, is 
compatible with the observed event yield in the signal region, Na = 138 events. 
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6 Systematic Uncertainties 



Table 4: Systematic uncertainties in the dilepton channels for given mass point and optimized 
mass window for the background (columns 3-6) and signal (columns 7-8) expected yields, 
following the procedures described in the text. In addition to the estimated signal uncertainties 
listed in the Table, constant uncertainties are considered on the integrated luminosity (2.2%), 
the lepton reconstruction and trigger efficiencies as determined by the "tag and probe" method 
(2%) [56] and the V mass selection as determined from a sample of boosted tt events (9%). 





Background 


Signal 


Mass Point 


Mass Window 


Stat. 


Fit 


Diff. w/ 


JES 


JES 


PDF 


(GeV) 


(GeV) 


(%) 


variations (%) 


MC (%) 


(%) 


(%) 


(%) 


W model, electron channel 


700 


640-760 


8 


3 


1 


4 


1 


4 


800 


755-845 


8 


21 


1 


5 


1 


5 


900 


855-945 


9 


21 


2 


7 


1 


5 


1000 


930-1070 


11 


17 


5 


2 


1 


6 


1100 


1020-1180 


12 


22 


4 


3 


1 


6 


1200 


1130-1270 


15 


26 


6 


7 


1 


6 


1300 


1220-1380 


17 


46 


6 


41 


1 


7 


1400 


1320-1480 


20 


64 


28 


14 


1 


8 


1500 


1390-1610 


23 


72 


31 


26 


1 


9 


RS mode 


, electron channel 


750 


690-810 


8 


14 


3 


2 


1 


4 


1000 


940-1060 


11 


17 


5 


7 


1 


6 


1250 


1180-1320 


16 


36 


3 


8 


1 


7 


1500 


1390-1610 


23 


72 


3 


26 


1 


9 


1750 


1540-1960 


31 


64 


10 


48 


1 


10 


2000 


1760-2240 


42 


42 


26 


110 


1 


11 


W model, muon channel 


700 


640-760 


7 


10 


4 


6 


1 


4 


800 


755-845 


8 


19 


6 


8 


1 


5 


900 


855-945 


9 


16 


2 


8 


1 


5 


1000 


930-1070 


10 


15 


8 


1 


1 


6 


1100 


1020-1180 


12 


20 


8 


1 


1 


6 


1200 


1130-1270 


14 


29 


21 


2 


1 


7 


1300 


1220-1380 


16 


43 


30 


6 


1 


7 


1400 


1320-1480 


19 


50 


68 


10 


1 


8 


1500 


1390-1610 


22 


46 


54 


24 


1 


9 


RS model, muon channe 




750 


690-810 


7 


14 


2 


7 


1 


4 


1000 


940-1060 


10 


15 


6 


7 


1 


6 


1250 


1180-1320 


15 


37 


4 


1 


1 


7 


1500 


1390-1610 


22 


46 


54 


24 


1 


9 


1750 


1540-1960 


30 


33 


31 


52 


1 


11 


2000 


1760-2240 


40 


64 


23 


130 


1 


11 
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Table 5: Systematic uncertainties in the E™ 1SS channel for the expected signal yields for the W' 
mass range M w E [700, 1500] GeV and graviton mass range M Gkk S [750,2000] GeV. 



Mass Point (GeV) 


PDF (%) 


JES (%) 


E™ 1SS (%) 


Total (%) 


W' model 


700 


4 


9 


9 


13 


800 


5 


8 


8 


12 


900 


5 


8 


8 


12 


1000 


5 


7 


7 


11 


1100 


6 


5 


6 


10 


1200 


6 


2 


4 


7 


1300 


7 


1 


3 


8 


1400 


8 


1 


3 


9 


1500 


8 


1 


3 


9 


RS model 


750 


4 


7 


7 


10 


1000 


4 


1 


3 


5 


1250 


4 


1 


3 


5 


1500 


4 


1 


3 


5 


1750 


4 


1 


3 


5 


2000 


4 


1 


3 


5 



6.2 Signal systematics 

There are several systematic uncertainties in the expected signal yields that are common across 
channels. These uncertainties are on the luminosity measurement, the JES effects on jets, the 
PDF, and the trigger and reconstruction efficiencies. A value of 2.2% was taken for the uncer- 
tainty in the measurement of the integrated luminosity |[57| . 

To determine the effect of JES uncertainty, we vary the jet px to pj ± Pj E g (pj,t]), where oj E g ( p x , rj ) 
is the total jet uncertainty, and apply the full analysis selection. The differences in the signal 
yields N +S i g and N_ S i g with respect to the nominal selection N^g are taken as the ±1(7 esti- 
mates for the JES systematic uncertainty. For W' and RS signals with the mass in the range 
[700, 2000] GeV in the dilepton channels and in the range [1250, 2000] GeV in the E™ ss channel, 
this systematic uncertainty is less than 1%. However, for resonance masses in the range [700, 
1200] GeV in the E^ 1SS channel this systematic uncertainty is found to be between 2 and 9%, 
owing to threshold effects. To estimate the systematic uncertainty associated with the choice of 
the PDF used for the simulated samples, three scenarios are considered: CTEQ6.6, MSTW2008 
and NNPDF2.0 [58]. The systematic uncertainty is set to half of the difference between the 
maximum and the minimum PDF values predicted for each mass point [59] . 

6.2.1 Dilepton channels 

To account for differences in trigger and reconstruction efficiencies between the Monte Carlo 
simulation and data, we determine scaling factors by using data control samples of Z — > pp 
and Z — > ee candidate events [60. 61]. We derive corrections for the muon (0.974 ± 0.001) and 
the electron (0.960 ± 0.004) channels and we apply them to the expected signal yields. These 
numbers assume that the efficiency does not vary with pj (Ej). However, we observe a small 
decrease (increase) in the efficiency in the asymptotic high-pj (high-Ej) region for muons (elec- 
trons) of about 2%. This small difference is used as the systematic uncertainty in the expected 
number of signal events for each mass point considered in this study. Finally, we assign a 
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9% systematic uncertainty on the V mass selection efficiency. This is determined by studying 
an independent sample of boosted tt — > Wb Wb events in which one of the W bosons decays 
leptonically and the other hadronically. 

6.2.2 E!f iss channel 

Propagating the jet energy scale effects to the calculation of E™ 1SS , and accounting for the an- 
ticorrelation between jets and E™ 1SS itself, we estimate a systematic effect of around 3% for all 
values of Mq studied, except for the lowest Mq = 750 GeV. In this case, because of threshold 
effects, the systematic effect is found to be around 7%. 

Summing in quadrature the uncertainties above, we arrive at a final 5% systematic uncertainty 
on the signal acceptance and efficiency except for Mq = 750 GeV, where a value of 10% is 
obtained for the final systematic uncertainty on the signal acceptance and efficiency. 



7 Results 

We do not observe any significant excess over the expected background. We employ the mod- 
ified frequentist CLg statistical method |62l l63j to search for exotic VZ resonances. For the 
dilepton channels we use a series of search windows corresponding to different mass hypothe- 
ses. Each mass window is optimized to give the best exclusion limit, a procedure which is 
also appropriate for establishing a new resonance discovery. The mass windows optimiza- 
tion has been carried out separately for the W and RS graviton hypotheses to account for 
differences in the width and efficiencies. For the E™ ss channel we perform a single counting 
experiment in the Mj > 900 GeV and My > 70 GeV region. We calculate 95% confidence level 
(CL) exclusion limits on the combined products of the cross section times the branching ratio 
cr(pp -> W') x H(W' -> WZ) and <j(pp -> Gkk) x B{G K k -> ZZ) for the three final states un- 
der study (separately and combined) as a function of the mass of the hypothetical resonance. 
We interpret these exclusion limits in two benchmark signal models: SSM W and RS graviton. 

The limit setting is performed by looking for an excess over the expected background in the VZ 
mass distributions for the three channels separately. Tables [6] and [7] show the search windows 
for each mass point with the corresponding signal efficiency and the numbers of expected back- 
ground and observed events in the electron and muon channels, respectively. These numbers 
are used as input for the calculation of the expected and observed exclusion limits on cross sec- 
tion times branching ratios at 95% CL that are also reported in the same Tables. Table [8] shows 
the signal efficiency and the observed and expected exclusion limits as a function of the signal 
mass in the E™ ss channel. The combined results are reported in Table [9 The exclusion limits 
as a function of the VZ resonance mass can be seen in Fig. |4j where a linear interpolation is 
used between the benchmark mass values. These limits can be interpreted in the theoretical 
framework of the W and RS graviton models. We exclude SSM W bosons with masses in the 
range 700-940 (890) GeV in the SSM at NNLO (LO) at 95% CL These results are complementary 
to the ones obtained in the tri-lepton analysis (with M w > > 1143 GeV in the SSM 1 64 ] ) . The ex- 
clusion limit calculated in the RS graviton model rules out masses (Mg kk ) in the range 750-880 
(800) GeV for k/ Mpi = 0.05 at NLO (LO). Assuming the resonance width is much smaller than 
the experimental resolution for the range of k/ Mpi considered here, the limit can be translated 
into the Mc KK -k/ Mpi plane. We do this by using the quadratic dependence of the cross section 
on k/ Mpi, and by assuming that the signal efficiency remains the same. The result is shown in 

These results are particularly relevant in the context of RS models proposed in recent studies 
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[65) , with SM fields propagating in the extra dimension where the graviton coupling to light 
fermions is strongly suppressed. This opens the possibility to an enhancement of the branching 
fractions for final states with V pairs, and motivates the investigation of large values of k/ Mpi 
II65I . In this scenario the previously-published searches for RS gravitons decaying to 77 and 
l + £~ final states [24] do not impose stringent bounds, since the branching fraction for these 
final states is suppressed. The results derived in this analysis are currently the most stringent 
in the V pair channel, and provide important constraints that are complementary to the ones 
from the search for resonances decaying to boosted top pairs ||66| . 

Table 6: Electron channel: Search window for each mass point with the corresponding sig- 
nal efficiency ("e S ig") and the numbers of mean expected background ("Nb g d") and observed 
("N b s ") events. The uncertainties include both statistical and systematic effects. These num- 
bers are used as input for the calculation of the expected and observed exclusion limits on 
tr(pp -> W') x H(W -> WZ) and a(pp -> G KK ) x B(G KK -> ZZ) at 95% CL which are re- 
ported in the last two columns. 



Mass Point 
(GeV) 


Window 
(GeV) 


N bgd 


N obs 


e sig 
(%) 


Obs. Limit 
(Pb) 


Exp. Limit 
(Pb) 






w' 


model 








700 


640-760 


39.7 ± 3.9 


43 


37±4 


0.44 


0.37 


800 


755-845 


24.6 ± 5.7 


23 


36±4 


0.33 


0.35 


900 


855-945 


17.1 ± 4.2 


12 


40±4 


0.18 


0.24 


1000 


930-1070 


17.1 ± 3.5 


17 


49±5 


0.20 


0.20 


1100 


1020-1180 


12.0 ± 3.0 


13 


48±5 


0.20 


0.18 


1200 


1130-1270 


6.3 ± 1.9 


5 


41±5 


0.13 


0.15 


1300 


1220-1380 


4.4 ± 2.8 


6 


32±4 


0.25 


0.20 


1400 


1320-1480 


2.7 ±1.8 


2 


23±3 


0.18 


0.19 


1500 


1390-1610 


2.5 ± 2.0 


2 


19±2 


0.22 


0.22 






RS 


model 








750 


690-810 


37.1 ± 6.0 


32 


27±3 


0.21 


0.25 


1000 


940-1060 


14.6 ± 3.1 


16 


35±4 


0.14 


0.13 


1250 


1180-1320 


4.9 ± 1.9 


7 


35±4 


0.11 


0.08 


1500 


1390-1610 


2.5 ± 2.0 


2 


27±3 


0.08 


0.08 


1750 


1540-1960 


2.0 ± 1.7 


1 


16±2 


0.10 


0.12 


2000 


1760-2240 


1.3 ± 1.6 





17±2 


0.06 


0.05 



8 Summary 

A search for new exotic particles decaying to the VZ final state was performed, where V is 
either a W or a Z decaying to hadrons, and the Z decays to electrons, muons, or a neutrino 
pair. The analysis is based on a data sample of pp collisions corresponding to an integrated 
luminosity of 5.0 fb _1 collected by the CMS experiment at the LHC at ^ = 7TeV in 2011. No 
significant excess is observed in the mass distribution of the VZ candidates compared with the 
background expectation from standard model processes. Lower bounds at the 95% confidence 
level are set in two theoretical models on the mass of hypothetical particles decaying to the 
VZ final state. Assuming heavy charged vector bosons in the sequential standard model, W 
bosons are excluded with masses in the range 700-940 (890) GeV at NNLO (LO). In the Randall- 
Sundrum model, graviton resonances with masses in the range 750-880 (800) GeV at NLO (LO) 
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Table 7: Muon channel: Search window for each mass point with the corresponding sig- 
nal efficiency ("e S ig") and the numbers of mean expected background ("Nb g d") and observed 
("N bs") events. The uncertainties include both statistical and systematic effects. These num- 
bers are used as input for the calculation of the expected and observed exclusion limits on 
cr(pp -> W') x H(W' -> WZ) and a(pp -> G KK ) x 23(G K k -> ZZ) at 95% CL, which are 
reported in the last two columns. 



Mass Point 
(GeV) 


Window 
(GeV) 


N bgd 


N bs 


6 sig 
(%) 


Obs. Limit 
(Pb) 


Exp. Limit 
(pb) 






W' 


model 








700 


640-760 


48.7 ± 8.9 


45 


40±4 


0.40 


0.45 


800 


755-845 


28.6 ± 6.9 


21 


40±4 


0.25 


0.34 


900 


855-945 


19.2 ± 4.3 


23 


41±4 


0.37 


0.29 


1000 


930-1070 


18.7 ± 3.7 


26 


51±6 


0.34 


0.22 


1100 


1020-1180 


12.9 ± 3.1 


12 


52±6 


0.17 


0.18 


1200 


1130-1270 


6.7 ± 2.2 


8 


44±5 


0.18 


0.15 


1300 


1220-1380 


4.6 ± 2.1 


4 


42±5 


0.13 


0.13 


1400 


1320-1480 


2.9 ± 2.0 


1 


39±5 


0.08 


0.11 


1500 


1390-1610 


2.6 ± 1.7 


2 


40±5 


0.11 


0.11 






RS 


model 








750 


690-810 


44.1 ± 9.2 


34 


30±3 


0.19 


0.26 


1000 


940-1060 


15.9 ± 3.4 


20 


39±4 


0.17 


0.13 


1250 


1180-1320 


5.2 ± 2.1 


6 


41±5 


0.08 


0.07 


1500 


1390-1610 


2.6 ± 1.7 


2 


44±6 


0.05 


0.05 


1750 


1540-1960 


2.1 ± 1.4 


2 


32±5 


0.06 


0.06 


2000 


1760-2240 


1.3 ± 1.9 


2 


42±6 


0.06 


0.05 
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Table 8: E™ 1SS channel: Expected and observed exclusion limits on cr(pp — > W') x H(W —> WZ) 
and 0"(pp — » Gkk) x £>(Gkk — > ZZ) at 95% CL for each mass point with the corresponding 
signal efficiency ("e s ig")- m tne Mr > 900 GeV region the expected background is B est = 153 ± 
29, including both statistical and systematic uncertainties, and the number of observed events 
is 138. These parameters are common for all mass points considered in this channel. 



Mass Point (GeV) 


£sig (%) 


Obs. Limit (pb) 


Exp. Limit (pb) 


W' model 


700 


0.2±0.1 


29 


33 


800 


0.9±0.1 


7.0 


8.2 


900 


8.0±0.5 


0.77 


0.90 


1000 


31±2 


0.19 


0.23 


1100 


49±2 


0.13 


0.15 


1200 


58±3 


0.11 


0.13 


1300 


64±3 


0.10 


0.11 


1400 


66±3 


0.09 


0.11 


1500 


69±3 


0.09 


0.11 


RS model 


750 


0.7±0.1 


4.1 


4.8 


1000 


25±2 


0.12 


0.14 


1250 


43±3 


0.07 


0.08 


1500 


54±3 


0.06 


0.07 


1750 


60±3 


0.05 


0.06 


2000 


63±3 


0.05 


0.06 



Table 9: Combined channels: Expected and observed exclusion limits on c(pp — > W') x 
B(W -> WZ) and a(pp -> Gkk) x B (Gkk ZZ) at 95% CL for the electron, muon, and 
£miss cnanne l s combined for each mass point and search window. 



Mass Point (GeV) 


Window 


Obs. Limit (pb) 


Exp. Limit (pb) 


W' model 


700 


640-760 


0.30 


0.25 


800 


755-845 


0.14 


0.21 


900 


855-945 


0.19 


0.18 


1000 


930-1070 


0.19 


0.15 


1100 


1020-1180 


0.11 


0.12 


1200 


1130-1270 


0.09 


0.09 


1300 


1220-1380 


0.09 


0.09 


1400 


1320-1480 


0.05 


0.08 


1500 


1390-1610 


0.07 


0.08 


RS model 


750 


690-810 


0.12 


0.16 


1000 


940-1060 


0.13 


0.09 


1250 


1180-1320 


0.07 


0.05 


1500 


1390-1610 


0.04 


0.04 


1750 


1540-1960 


0.05 


0.05 


2000 


1760-2240 


0.04 


0.03 
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CMS Vs = 7TeV 



.Q 
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Q. 
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1 1 1 I 1 1 1 1 I 1 1 
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M wz (GeV) 
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N 

T 
o 

CO 



T 

Q. 
Q. 



I 1 1 1 I 1 1 1 I 1 1 1 I 
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Exp. + 1o [] Exp. ± 2a 

a RS (LO) - a RS (NLO) 




Limit RS (LO) = 803 GeV _ 
Limit RS (NLO) = 879 GeV 



800 1000 1200 1400 



1600 1800 2000 

M zz (GeV) 



Figure 4: Observed and expected 95% CL upper cross section limits and comparison with the 
theoretical predictions in W (top) and RS graviton with k/ Mj>\ = 0.05 (bottom) models for the 
combination of electron, muon, and EJp lss channels. The limits are calculated with the modified 
frequentist CLg statistical method. 
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CMS ^s = 7TeV 



t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

; ATLAS (ZZ 2l2j/4j, LO, 1.02/fb) 



Ldt = 5.0 fb 



CMS (ZZ 2l2j, LO, 4.9/fb) 
CMS (ZZ ^ 1j2l/1j2v, LO, 5.0/fb) 



CMS (ZZ -> 1j2l/1j2v, NLO, 5.0/fb) 

Precision EW observables 



0.3 



0.25 



0.2 



'0.15 



0.1 



0.05 



600 800 1000 1200 1400 

Graviton Mass (GeV) 

Figure 5: Observed exclusion limits from this analysis (ZZ — > Ij2£/lj2v) interpreted in the 
context of the RS graviton model, assuming a LO prediction (solid green line) and NLO pre- 
diction (dot-dashed green line). Also shown are results from ATLAS (ZZ — > 2£2j/Aj) II35I 
(dot-dot-dashed red line) and another CMS publication (ZZ — > 212]) 1)361 (dotted blue line). 
The shaded region corresponds to the indirect limits derived from precision electroweak ob- 
servables in 




are excluded for k/M?\ =0.05. These are the first results from the LHC on VZ searches using 
final states with a boosted massive jet and a lepton pair or missing transverse energy. 
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